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ABSTRACT 


This  report  describes  a  methodology  put  into  a  preprocessor 
program  for  computer  war-game  models.  The  preprocessor  computes 
kill  rates  or  probability  of  kills  (Pk)  for  selected  direct-fire 
weapons  against  selected  targets.  These  are  input  data  used  by 
the  war-game  models. 

The  methodology  is  an  expansion  to  a  methodology  already  in 
the  preprocessor.  This  methodology  computes  probability  of  hit 
(Ph)  for  the  weapon  fired  against  the  target.  Ph  is  an  inter¬ 
mediate  value  used  in  computing  the  kill  rate  or  Pk. 

The  preprocessor  needed  the  expansion  to  acco.-nmodate  a  new 
type  of  weapon  which  has  a  flight  path  above  the  target.  The 
previous  methodology  assumed  a  ground-level  flight  path. 


V 


1.  Introduction.  A  White  Sands  Missile  Range  (WSMR)  program 

called  IMrect-fi :”e  Weapon  Preprocessor  (DFWP)  computes  kill  rates 
for  Vector- In-Commandar  (VIC)  and  other  computer  war-game  models. 
The  kill  r:  data  are  for  selected  direct-fire  weapons  against 

selected  targets.  Because  of  new  top-attack  weapons,  the  method¬ 
ology  in  this  progj’am  was  not  adequate.  This  report  discusses 
the  problem  and  defines  the  solution. 

2.  Problem.  The  methodology  for  probability  of  hit  (Ph)  needed 
to  be  changed. 

a.  Overview  of  Ph  computation. 

(1)  DFWP  uses  two  kinds  of  data. 

(a)  One  is  the  target  dimensions.  This  is  available  from 
the  Ballistic  Research  Laboratory  (BRL) . 

(b)  The  other  is  the  probability  density  data  for  round 
impacts  of  the  weapon.  This  is  available  from  the  Army  Materiel 
Systems  Analysis  Activity  (AMSAA) . 

(2)  The  Ph  methodology  has  two  parts. 

(a)  DFWP  computes  target  area  from  the  target  dimensions.. 
and  weapon  flight  path.  The  orientation  of  the  target  (aspect 
angle)  and  vertical  angle  of  weapon  approach  (attack  angle) 
determine  the  flight  path. 

(b)  DFWP  computes  Ph  by  integrating  the  probability 
density  function  (pdf)  of  round  impact  p(X,Y)  over  the  target 
ar  ea , 


I  I 

Ph  =  !  I  p(X,Y)  dY  dX 

I  I 

Target  area 


b.  Problem. 

(1)  DFWP  must  compute  target  area  differently  for 
top-attack  weapons  than  for  the  other  weapons.  Those  weapons 
have  an  attack  angle  of  zero  degrees  to  the  ground.  The  target 
area  computes  to  be  a  rectangle  or  rectangles.  Target  area  for 
top-attack  weapons  must  be  cc  iputed  by  an  expanded  methodology, 
and  it  computes  to  be  an  irregular  polygon. 
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(2)  DFWP  must  compute  Ph  differently  for  top-attack 
weapons.  DFWP  computes  PH  for  the  other  weapons  (paragraph 
2a(2)(b))  by  table  lookup.  This  is  because  the  target  area  is  a 
rectangle  or  two  rectangles.  DFWP  must  compute  Ph  for  top-attack 
weapons  by  a  numerical  method  since  the  target  area  is  more 
oompl i Gated . 

3 .  Assumptions . 

a.  The  DFWP  assumes  the  targets  to  be  a  box  or  two  boxes, 
one  on  top  of  the  other.  If  two  boxes  are  assumed,  then  the 
bottom  box  is  called  the  hull,  and  the  top  box  is  called  the 
turret.  The  turret  is  assumed  to  be  exactly  on  the  center  of  the 
hull  (see  figure  i  (a)  and  (c)).  (This  last  assumption  might 
change  in  the  future.) 

b.  The  DFWP  assumes  p(X,Y)  to  be  normal  and  independent  in 
the  horizontal  and  vertical  directions. 

c.  The  DFWP  computes  both  fully  exposed  and  defilade 
targets.  Defilade  targets  are  behind  a  barrier  of  sorts  and  are 
partially  hidden  from  view.  The  distance  of  the  barrier  to  the 
target  in  addition  to  the  barrier  height  determines  the  exposed 
part  of  the  target  for  top-attack  weapons.  Barrier  height  is  an 
input  to  DFWP  but  not  barrier  distance.  DFWP  as.sumes  the  barrier 
distance  equal  to  the  hull  height. 

d.  The  DFWP  assumes  aimpoint  to  be  the  center  of  target 
volume . 
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4 .  Target  area . 


a .  General . 

(1)  Summary.  The  DFWP  computes  target  area  from  the 
weapon  flight  path  and,  the  target  dimensions.  The  aspect  angle 
and  attack  angle  determine  the  weapon  flight  path.  DFWP  uses  a 
coordinate  system  for  the  target  dimensions. 

(2)  Aspect  angle.  The  DFWP  computes  Ph  for  the  weapon 
engaging  the  target  from  the  front,  side,  and  other  directions. 
This  direction  is  the  aspect  angle.  Engagement  from  the  front  is 
a  cero-degree  aspect  angle.  Engagement  from  the  side  is  a  90- 
degree  aspect  angle. 

(3)  Attack  angle.  Each  weapon  has  an  attack  angle. 
Paragraph  b  discusses  weapons  with  an  attack  angle  of  zero 
degress.  This  was  the  only  case  considered  previously,  and 
paragraph  b  discusses  the  methodology  that  was  already  in  the 
DFWP.  Paragraph  c  discusse.s  weapons  with  an  attack  angle  of  90 
degrees.  Paragraph  d  discusses  weapons  with  an  attack  angle 
between  0  and  90  degrees. 

^4)  Coordinate  system.  The  corner  coordin-ates  describe 


the  target  and  barrier. 

The 

corners  are: 

(See  f igure  1  (b) ,  (d) , 

and 

(f )  .  ) 

Hull 

Turret 

Barrier 

1  -  Front  bottom 

2  -  Back  top 

6 

Back  top 

10  ■*  Front  top 

3  -  Left  top 

7 

Left  top 

4  -  Right  top 

8 

Right  top 

The  coordinate  system  has  the  X-axis  for  the  horizontal 
direction,  the  Y-axis  for  the  vertical  direction,  the  Z-axis  for 
the  depth  direction,  and  the  origin,  0,  at  the  center  top  of  the 
hull.  The  target  is  rotated  about  the  origin  for  the  different 
aspect  angles.  The  coordinates  for  point  1  are  XI,  Y1 ,  Zl,  and 
W1  ;  the  coordinates  for  point  2  are  X2 ,  Y2 ,  Z2 ,  and  W2 ;  etc.  The 
X  and  Y  coordinates  describe  the  target  area  for  level  attack. 
The  X  and  Z  coordinates  describe  the  target  area  for  straight 
down  attack.  X  and  W  coordinates  describe  the  target  area  for 
angle  attack.  The  W-axis  is  an  additional  axis  made  perpen¬ 
dicular  to  the  round  flight  path  and  to  the  X-axis. 

(5)  Figure  1  shows  target  areas.  The  dimensions  are: 

W  =  Hull  width  =  2  meters  (m)  w  =  Turret  width  =  .5m 

L  =  Hull  length  =  3  m  1  =  Turret  length  =  1.5  m 

H  =  Hull  height  =  1.75  m  h  =  Turret  height  =  2  m 

Hb  =  Barrier  height  =  1.5  m 
A  =  Aspect  angle  =  30  degrees 
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LeV'-l  attack. 


( ’  i  j  Target  area  formulas. 

(a)  Coordinates  for  hull  corners  3  and  4. 

X4  =  W  /2  »  cos(A)  +  L  /2  »  sin(A) 

X3  "  -X4 
y3  =  y4  =  0 

(b)  Coordinates  for  turret  corners  7  and  8. 

X8  =  w  -■  2  *  cos  (A)  +  1  /2  *  sintA) 

X7  =  -XS 
y7  Y8  =  h 

(c)  Coordinates  for  barrier  corner  10. 

XIO  =  -  (W  /2  +  H)  if  cos(A)  +  (L  /2  +  H)  *  sin(A) 
YIO  =  -H  +  Hh 


(d)  Coordinates  for  aimpoint. 

XO  =  0 

,  YO  =  (h  +  YIO)  /  2 

(2)  Ph  formula.  DFWP  integrates  pdf  p(X,Y)  over  the 
target  area  of  hull  and  turret. 


Ph  = 

X4 

1 

Ph  on 

0 

) 

hul  1  + 

Ph  on 

X8  h 

1  1 

turret 

Ph  =  ! 

1  p(X.Y) 

dY  dX  + 

'  '  p(X,Y) 

dY  dX 

I  I  I 


X3  YIO  X7  0 
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straight-down  attack. 

(1)  Target  area  formulas. 


0  . 


(a)  Hull  coordinates. 


Let:  A1  =  A  =  aspect  angle. 

A2  =  A1  +  90. 

M  =  rotation  matric. 

M  '  jcos(Al)  -sin(Al) 
isin(Al)  cos^Al) 

I  '  I 

!  L/  2  1 

I  W/2| 
I-L/21 


Then 


XI 

=  M 

-W/2  i 

j  X2 

j  =  M 

Z1 

-L/2! 

1  Z2 

X3 

=  M 

-W/2| 

j  X4 

=  M 

Z3 

L/2I 

1  24 

(b)  Turret  coordinates.  The  turret  area  is  not  computed. 
The  turret  area  is  inside  the  hull  area. 


fc)  Earner  coordinates.  The  barrier  is  not  computed. 
The  barrier  does  not  cover  the  target  area. 


(d)  Aimpoint  coordinates. 

XO  =  I  0 
ZO  I  9 


(e)  Equations  for  the  sides  (See  figure  2) . 


L  1  f  X ) 

=  ZI 

4 

tan ( A2 ) 

(X  - 

XI ) 

L2  (X) 

=  Z2 

4> 

tan ( A ] ) 

(X  - 

X2) 

IM  (X) 

=  7,1 

-f 

tan ( A1 ) 

# 

(X  - 

XI  ) 

H  2  (X  /' 

—  '■70 

—  w  ^ 

^  f  A  n  > 
v/ci-ii  \  nxf  / 

J 

K 

<  AT 

X2) 

6 


zt 


L2  (X) 


I 


Figure  2.  Straight-down  area 


(2)  Ph  formula. 

If  XI  <  X2  (Shown  in  figure  2.) , 


(Section  1) 

XI  L2(X) 

I  I 

Ph  =  !  !p(X,Z)  dZ  dX 

I  I 

X3  Ll(X) 

If  X2  <  XI, 

( Sect i on  1 ) 

X2  L2(X) 

I  I 

Ph  =  I  ip(X,Z)  dZ  dX 

I  I 

X3  L1(X) 


(Section  2) 

X2  L2(X) 

I  I 

+  I  'p(X.Z)  dZ  dX 

I  j 

XI  RKX) 


(Section  2) 

XI  R2(X) 

I  I 

+  !  'p(X.Z)  dZ  dX 

I  I 

X2  LUX) 


(Section  3) 
X4  R2(X) 


+ 


p(X.Z)  dZ  dX 


X2  RKX) 


(Section  3 ) 

X4  E2(X) 

I  I 

+  '  ip(X.Z)  dZ  dX 

I  I 

XI  RKX) 


(c  ) 


Weapon  vien 

Figure  3.  Angle  area 


8 


d.  Angle  attack, 


(1) 

(a) 

1. 

a . 


b  . 


c  . 


d  . 


2  . 

a  . 

b  . 
3  . 
a  . 
b  . 


Target  area  formulas. 

Coordinates . 

DFWP  first  computes  the  X-2  plane. 


Let:  A1  =  aspect  angle 

A2  =  A1  +  00 
M  =  rotation  matrix 


cos(Al)  -sin(Al) 
sin(Al)  cos(Al) 


Hull  corner  coordinates. 


XI 

-  W/  2  j 

1X2 

=  M 

W/2 

Z1 

=  M  j 

-L/2  1 

Z2 

L/2 

X3 

=  M 

-W/2  1 

X4  1 

=  M 

W/2 

i 

23 

L/2  1 

Z4  1 

!  -L/2 

Tur  r 

'et 

corner 

coordinates. 

X6  j 

w/2 

26  1 

j  =  M 

1/2 

X7| 

=  M 

-w/2j  1 

X8  1 

=  M 

w/  2 

27! 

1/2  1 

23 

-1/2 

Barrier  corner  coordinates. 


XIO 

=  M 

- (W/2+H) 

2  10 

-  (L/2^-H) 

DFWP  then  computes  the  Y-Z  plane, 

Y  coordinates.  See  paragraph  4b(l/ . 

2  coordinates.  See  paragraph  4d  (  1 )  (a)  1_. 

DFWP  then  computes  the  X-W  plane  (target  area  plane) , 


X  coordinates.  See  paragraph  4d ( 1 ) (a)^. 

W  coordinates.  Figure  3  (b)  shows  the  side  view  of 
the  target  (Y-Z  plane).  It  shows  the  angle  ,Ai  of 
flight  path  and  the  W-axis  which  is  the  edge  of  X-W 
plane  perpendicular  to  the  flight  path.  The  Y  and  Z 
cooi'dinates  and  A  compute  the  W  coordinate  by: 


W7.  =  Z*/.  sin  (A)  +  Y7,  *  cos  (A). 
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(b)  Barrier  coverage.  The  methodology  computes  the  hull 
area  that  the  barrier  does  not  cover. 

1_.  Ways  barrier  covers  hull.  The  position  of  the 

barrier  corner  10  in  relation  to  the  bottom  edges  of 
the  hull  determines  the  way  the  barrier  covers  the 
hul  1  . 

-  Case  1.  Point  10  is  below  both  edges.  Neither  side 
of  barrier  covers  hull. 

-  Case  2.  Point  10  is  above  both  edges.  Both  sides  o 
barrier  cover  hull. 

-  Case  3.  Point  10  is  above  right  edge  and  below  left 
edge.  Right  side  of  barrier  covers  hull.  (This  is 
the  case  in  figure  3  (c)  .) 

-  Case  4.  Point  10  is  below  right  edge  and  above  left 
edge.  Left  side  of  barrier  covers  hull. 

2.  Determination  of  way  barrier  covers  hull.  Points  1 
and  10  get  U  and  V  coordinates.  The  U-saxis  is 
perpendicular  to  the  left,  bottom  edge  of  hull,  and 
the  V-axis  is  perpendicular  to  the  right  bottom  edge 
of  hull.  See  figure  4. 


A 


I 


I 


Figure  4.  U  and  V  axes 
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Let  U  be  angle  of  U  axis  and  V  be  angle  of  V  axis. 


Then  : 

cos ( U) 
s in (U) 
cos (V) 
sin ( V) 

=  cos (B2-90 )  = 

=  sin(B2-90}  =  ■ 

=  cos ( B  1  +  90 )  =  ■ 

=  sin(Bl+90)  = 

sin (B2) 
-cos ( B2 ) 
-sin (B1 ) 
cos (B1 ) 

Let : 

M 

=  1  cos (C)  s in (U)  I 
|cos(V}  sin(V)l 

Then  ; 

juiol  = 

(Viol 

M  1 XlO  j  1 U1 j  = 

1 WIO 1  ! VI 1 

M  XI 

W1 

-  Case  1  . 

UlO  < 

U1  and  VIO  <  VI . 

-  Case  2 . 

UlO  > 

U1  and  VIO  >  VI. 

-  Case  3  , 

UlO  < 

U1  and  VlO  >  VI. 

-  Case  4. 

UlO  > 

U1  and  VIO  <  VI. 

3.  Computation.  Point  1,  the  bottom 
changed  to  the  bottom  point  of  the 
the  barrier  does  not  cover. 

point  of 
part  of 

hull,  13 
hull  th.a.t 

"  Case  1. 

Point 

1  is  not  changed. 

-  Case  2. 

Point 

1  is  changed  to  the 

barrier 

point  10. 

-  Case  3. 
i  s  the 
bottom 

Point  1  is  changed  to  the  point  11.  Point 
intersection  of  right  top  edge  of  barrier  with 
left  edge  of  hull.  Coordinates  of  point  11 

are  ; 

XI  1  = 

(tan(B2)»  XI 

-  tan ( B  1 )  + 

XIO 

+  WIO  - 

W1  ) 

/  (tan(B2)  - 

tan(Bl) ) 

W1 1  = 

(cot (B2) «  W1 

-  cot (B  1 ) » 

WIO 

+  XIO  - 

XI ) 

/  ( c  0 1  ( B  2 )  - 

cot  (Bl)  ) 

Case  4. 

Point  1  is 

changed  to 

the 

point  12 

Point 

12  is  the  intersection  of  left  top  edge  of  barrier  with 

bottom  right  edge  of  hull.  Coordinates  of  point  12 
are  : 

X12  =  (tan(B2)»  XIO  -  tan(Bl)»  XI  +  W1  -  WIO) 

/  (tan(B2)  -  tan(Bl)) 

W12  =  (cot(B2)»  WIO  -  cot(Bl)»  W1  +  XI  -  XIO) 

/  (cot  (B2)  -  cot (B1 ) ) 


1  1 


(c)  Turret  coverage.  The  methodology  computes  the  turret 
area  that  does  not  cover  the  hull. 


area 


1.  Ways  turret  extends  over  hull.  The  position  of  point 
6  in  relation  to  the  top  edges  of  the  hull  determines 
the  way  the  turret  extends  over  the  hull. 


-  Case  1 .  Po i n t 

6 

is  below 

both 

edges  . 

The  turret 

is  within  the 

hull  area. 

-  Case  2.  Point 

6 

is  above 

both 

edges  . 

The  turret  may 

extend  over  both 

edges  . 

-  Case  3 .  Po i n t 

6 

is  below 

right 

edge 

and  above 

left 

edge.  The  turret 

extends 

above 

left 

edge  . 

~  Case  4 .  Point 

6 

is  above 

right 

edge 

and  below 

left 

edge .  The  turret 

extends 

above 

right 

edge  . 

2_.  Determination 

o  f 

way  turret  extends  over  hull. 

Point 

6  and  point  2 

get 

U  and  V 

coordinates 

Let  U. 

V,  and 

M  be  the  same 

as 

in  paragraph  4d ( 1 ) (b)^. 

Then  ;  j  i 

S 

M  j  Xd  j 

iU2  j 

=■  M 

1X21 

1  V6  1 

1  wel 

1  V2  1 

1  W2l 

-  Case  1 .  U6  < 

U2 

and  V6  < 

V2  . 

-  Case  2.  U6  > 

U2 

and  VS  > 

V2. 

-  Case  3 .  U6  < 

U2 

and  VS  > 

V2  . 

-  Cas e  4 .  U6  > 

U2 

and  VS  < 

V2  . 

•3  .  Co  nip  Li  w  21  tL  i  o  n 

T-U  ^ 
a  11 

X  C  0  O  i' u  i  n  <si  0  e 

of  a  point  D,  wnich  is 

on  the  bottom  part  of  turret  not  covering  hull,  is 
computed . 


Case  1.  X5  is  not  computed. 


-  Case  2.  If  X7  <  X2  <  X8  (Turret  extends  over  both 
edges  of  hull.),  then  X5  is  made  equal  to  X2 .  If  X8 
<  X2  (Turret  extends  over  left  edge  of  hull.)  ,  then 
X5  is  made  equal  to  X8 .  If  X2  <  X7  (Turret  extends 
over  right  edge  of  hull.) ,  then  X5  is  made  equal  to 
X7  , 


12 


Case  3.  X2 1  is  computed.  Point  21  is  the 
intersection  of  top  right  edge  of  turret  with  top 
lef  t  edge  of  hull . 

X21  =  (tan(B2)»  X2  -  tan(Bl)»  X6  +  WS  -  W2 ) 

/  (tan(B2)  -  tan(Bl)) 

If  X8  <  X2 1  (Top  edge  of  turret  does  not  intersect 
with  top  edge  of  hull.)  ,  then  X5  is  mads  equal  to  X8 . 
Otherwise  X5  is  made  equal  to  X21. 

Case  4.  X22  is  computed.  Point  22  is  the 

intersection  of  top  left  edge  of  turret  with  top 
right  edge  of  hull. 

X22  =  (tan(B2)»  X6  -  tan(Bl)*  X2  +  W2  -  W6 ) 

/  (tan (B2)  -  tan (Bl) ) 

If  X22  <  X7  (Top  edge  of  turret  does  not  intersect 
with  top  edge  of  hull.),  then  X5  is  made  equal  to  X7 . 
Otherwise  X5  is  made  equal  to  X22. 

Angles.  A1  and  A2  are  labeled  Bl  and  B2  in  the  X-W 
plane.  See  figure  3  (c) .  Figure  5  shows  the  cosines 
and  sines  of  A1  and  A2  in  the  X-Z  plan©  and  of  El .and 
B2  in  the  X-W  plane.  Z  distances  are  changed  to  W 
distances  by  the  formula 

W  -  Z  *  sin(X) 

Therefore  the  formulas 

Dl  =  sqrt  (cos(Al)»f*2  +  (  s  i  n  ( A )  <^s  i  n  ( A 1  )  2  ) 
cos(Bl)  =  cos(Al)  /  Dl 

sin(Bl)  =  sin(A)  »  sin(Al)  /  Dl 
tan(Bl)  =  sin(Bl)  /  cos(Bl) 

D2  =  sqrt  (aos(A2)»^2  +  ( s i n ( A) *3 i n ( A2 ) « * 2 ) 
cos(D2)  =  cos(A2)  /  D2 

sin(B2)  =  sin(A)  »  sin(A2)  /  D2 
tan(B2)  =  sin(B2)  /  cos(B2) 

compute  Bl  and  B2. 
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"T 

s in ( A1 ) 

Y 


Figure  5.  Angles  of  target  sides 


Equations 

0  f 

the  sides. 

(See 

figure 

LI  (X) 

5 

W1 

+ 

tan (B2) 

tt 

(X  - 

XI ) 

L2  (X) 

W2 

-t- 

tan ( B 1 ) 

(X  - 

X2) 

L3  (X) 

W6 

+ 

tan  (B1 ) 

(X  - 

xe) 

Ri  (XJ 

W1 

+ 

tan  (B 1 ) 

* 

(X  - 

XI) 

R2(X) 

W2 

+ 

tan (B2) 

if 

(X  - 

X2) 

R3  (X) 

W6 

+ 

tan(B2) 

n 

(X  - 

xe) 
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(2)  Probability  of  hit  formulas. 

(a)  Ph  on  hul 1 . 

If  XI  <  X2  (This  is  the  case  in  figure  3.) 


(Section  1) 
XI  L2(X) 


(Section  2) 
X2  L2(X) 

I  I 


f(X)  !  f(Z)  dZ  dX  +  !  f(X)  !  f(2)  dZ  dX 


LI  (X) 

(Section  3) 


R1  (X) 


X4 

1 

R2  (X) 

1 

+  !  f(X)  1  f (Z) 

1  1 

X2 

El  (X) 

If 

X2  <  XI 

X2 

1 

L2  (X) 

( 

Ph 

=  !  f(X) 

1  f(Z)  dZ 

LL(X) 

X4 


(X)  XI  R2(X) 

I  ! 

f (Z)  dZ  dX  +  I  f (X)  i  f (2)  dZ  dX 

I  I 

X2  LUX) 

B2(X) 


+  I  f (X)  !  f (Z)  dZ  dX 

I  ! 

XI  RKX) 


(b)  Ph  on  turret. 

If  X5  <  X6. 

(Section  1) 
X5  L3(X) 


L2  (X) 


(Section  2) 
X6  L3(X) 


f(2)  dZ  dX  +  i  f(X) 


f (Z)  dZ  dX 


K2  (X) 


(Section  3) 
X8  R3(X) 


f (Z)  dZ  dX 


R2  (X) 
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If  X6 


Ph 


X5 


(This  is  the  case  in  figure  3.)  , 


X6 


L3  (X) 


X7 


!  f (X)  !  f (Z)  dZ  dX 

! 

L2  (X) 

X8  B3(X) 

! 

+  !  f(X) 

I 

X5  R2(X) 


X5  R3(X) 

I  ! 

!  f (X)  !  f (Z)  dZ  dX 

I  I 

X6  L2(X) 


f (Z)  dZ  dX 


(c)  Ph  on  target . 

Ph  =  Ph  on  hull  +  Ph  on  turret 


5 .  Probability  of  hit. 


a.  Level  attack.  Each  double  integral  in  the  Ph  formula 
from  paragraph  4b(2)  is  able  to  be  changed  into  a  product  of  two 
single  integrals,  thus; 


X4  0 

I  I 

Ph  =  I  f(X)  dX  *  I  f(Y)  dY 

!  I 

X3  YIO 


X8 

i 

+  I  f(X)  dX 
! 

X7 


h 

I 

»  I  f(Y)  dY 
! 

0 


Each  integral  can  be  evaluated  from  two  lookups  in  the  normal 
distribution  table. 


b.  Angle  or  s tr a i gh t -do wn  attack.  From  paragraphs  4c  (2)  and 
4d(2) ,  the  Ph  formula  for  a  section  is: 


YV 


Ph 


!  f  (X) 


f (Y)  dY  dX 


X7. 


K1  (X) 


The  outer  integral  must  be  evaluated  numerical ly .  The  inner 
integral,  though,  for  an  X  value  can  be  evaluated  from  two 
lookups  in  the  normal  distribution  table.  The  numerical 
integration  method  that  DFWP  implements  is:  First  DFWP  makes  a 
change  of  variable,  dU  =  f(X)  dX.  This  gets  rid  of  f(X)  in  the 
integrand.  The  integral  becomes: 
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UX  K2(X) 

I  I 

Ph  =  I  I  f(Y)  dY  dU 

(  I 

UX  KKX) 

The  normal  distribution  table  gives  the  outer  integral  limits. 

X*/. 

! 

U%  =  !  f(T)  dT  =  P(X%) 

I 

-  oo 

The  inverse  normal  table  gives  the  Xs  of  the  inner  integral 
limits . 


X  =  F'"  (U) 


The  DWFP  then  does  the  numerical  evaluation  using  Us  with  values 
.005,  .015,  .025,  etc.,  .995,  which  are  between  the  outer 

integral  limits.  These  values  are  midpoints  of  intervals  of 
width  .01  which  are  [0.,.01],  [.01, .02],  [.02, .03],  etc,, 

[.99,1.].  The  beginning  interval  and  last  interval  of  the 
numerical  integrati-on  will  not  necessarily  have  widths  of  .01. 


a . 

Input  values , 

W  = 

L  = 
H  = 
Hb 

Hull  width  =  2  m 

Hull  length  =  3  m 

Hull  height  =  1.75  m 
=  Barrier  height  =  1.5  m 

w  =  Turret  width  =  .5  m 
)  =  Turret  length  =  1.5  m 
h  =  Turret  height  =  2  m 

A1 

=  Aspect  angle 

= 

30 

A2  =  A1  +  90  =  120 

XE 

XD 

=  X-bias  =  0.2 
=  X-dispersion 

m. 

2  m 

WB  =  W-bias  =  0. 

WD  =  W-dispersion  =  1  m 

b  . 

Level  attack. 

(1)  Target  area. 

(a)  Hull. 

X4  =  W 
=  1  . 

/2 

* 

*  cos ( A 1 ) 
cos (30)  + 

+  L  /2  sin(Al) 

1.5  *  sin(30) 

X3 


1,616  m 
-X4 

-1,616  m 
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(b)  Turret, 

Xa  =  w  /2  cos(Al)  1  /2  «  sin(Al) 

=  0.25  *  cos(30)  +  0.75  it  sin(30) 

=  0.5915  m 

X7  =  -X8 

=  -0 .5915  m 

(c)  Barrier. 

YIO  =  -H  +•  Hb  =  -0.25  m 

(d)  Aimpoint. 

XO  =  0 . 

YO  =  (h  +  YIO)  /  2 

=  0 . 875  rn 

(  2 )  Ph  I  n  tej^ral  , 

h 

I 

*  !  f(Y)  dY 

t 

0 


X4  0  X8 

'  '  ! 

Ph  =  I  i(X)  dX  »  I  f(Y)  dY  +  '  f(X)  dX 

!  I  i 

X3  YIO  X7 


The  integral  limits  must  be  changed  to  standard  deviations  by  the 
formulas:  X-deviation  =  (X  -  XM)  /  XD  and  the  same  for  Y. 


The  mean  is  the  sum  of  aimpoint  and  bias. 
XM  =  XO+XB  =  0.+0.2  =  0.2m 

YM  =  YO  +  YB  =  .375  ■(-  0  .  =  0.875  m 


First  upper  limit  -  (X4  -0.2)  /  2.  =  0.70S  sd 
First  lower  limit  =  (X3  -0.2)  /  2.  =  -0.908  sd 


Second  upper  limit 
Second  1-ower  limit 


(0  -0.875)  /I.  =  -0.875  sd 

(YIO  -0.875)  /I.  =  -1.125  sd 


Third  upper  limit 
Third  lower  limit 


(X8  -0 . 2)  /  2 .  =  0  .  19575  Sd 
(X7  -0.2)  /  2.  =  -0.39575  Sd 


Fourth  upper  limit 
Fourth  lower  limit 


(h  -0.875)  /  1 .  =  1 . 125  sd 
(0  -0.875)  /I.  =  -0.875  sd 


Using  the  normal  distribution  table  F(X). 

Ph  below 
hull 

bottom  Side 
F(-I  .  125)  ) 


Ph 


Ph  left 
of  hull 
right  side 
(  F(0.708)  - 


Ph  left 
of  hull 
left  side 
F ( -0 . 908) ) 


Ph  below 
hull 

top  side 
(F(-0. 875) 


13 


Ph  left  Ph  left  Ph  below  Ph  below 

of  turret  of  turret  turret  turret 

right  side  left  side  top  side  bottom  side 

+  (F(0. 19575)  -  F(-0. 39575))  »  (F(1.125)  -  F(-0.875j) 

Ph  between  hull  Ph  between  hull 

right  and  left  sides  top  and  bottom  sides 

Ph  =  (.76053  -  .18194)  »  (.19079  ~  .13029) 

Ph  between  turret  Ph  between  turret 

right  and  left  sides  top  and  bottom  sides 

+  (.57760  -  .34614)  »  (.86971  -  .19079) 

Ph  =  .57859  »  .06050  +  .23146  »  .67892 

=  .19215 

c.  Straight-down  attack. 

(1)  Target  area. 

(a)  Coordinates . 


In 

meters . 

M 

=  1 

gos(30)  -si 

.n (30) 1 

1 

sin(30)  cos(30)i 

XI 

= 

M  -W/2j 

=  cos30»  (  -  1  .  ) 

sin30» ( • 

1.5) 

21 

-L/2I 

sin30»  (-1  . ) 

+ 

cos30» ( - 

1.5) 

= 

-0 .1160  m 
-1 . 7990  m 

X2  1 

M  W/2 

=  cos30  *  1  . 

sin30  * 

1  .  5  1 

22  1 

L/2 

s  in30  1  . 

+ 

cos30  * 

1.51 

= 

1  ^  1  ^  -  1 

V  ,  i  1 u  u  mi 

1  1 . 7990  m| 

X3 

I 

M  -W/2 

=  cos30»  (  -  1  .  ) 

sin30  * 

1  .  5 

23 

1 

1 

L/2 

s in30*  (  -  1  .  ) 

+ 

cos30  * 

1 . 5 

= 

-1.3160  m 

0.7990  m 

X4 

= 

M  W/2  1 

=  COS30  *  1  . 

_ 

s  in30» ( - 

1 .5) 

24 

-L/2  1 

s in30  *  1  . 

cos30»  ( ~ 

1  .5) 

1 . 6160  m 
-0 . 7990  m 
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In  standard  deviations 


(b)  Equations  of  sides 
1  .  In  meters . 


XI 

= 

(-0. 116  -XM) 

/XD 

(--0,116 

-0. 

Z1 

(-1.799  -ZM) 

/ZD 

(-1.799 

-c 

-0,158  sd 
-1.799  sd 

X2 

(0.116  -XM) 

/XD 

— 

-0 .042 

sd 

Z2 

(1.799  -ZM) 

/ZD 

1 . 799 

sd 

X3 

(-1.616  -XM) 

/XD 

=: 

-0 . 908 

sd 

Z3 

(0.799  -ZM) 

/ZD 

0.799 

sd 

X4 

(1.616  -XM) 

/XD 

Si 

0 . 703 

sd 

Z4 

(-0.799  -XM) 

/ZD 

-0 . 799 

sd 

/  1 


L1(X)  =  Z1  +  tan(A2)  *  (X-Xl) 

L2(X)  =  Z2  +  tan(Al)  #  (X-X2) 

RKX)  =  Zi  +  tan(Al)  «  (X-Xl) 

R2(X)  =  Z2  +  tan(A2)  »  (X-X2) 


2_.  In  standard  deviations.  The  cross  section  area  is 
different  when  the  X  and  Z  are  measured  in  standard 
deviations.  The  angles  A1  and  A2  are  different,  and, 
therefore,  tan(Al)  and  tan(A2)  are  different.  It  can 
be  shown  that  tan(Ai»)  =  XD  /  ZD  »  tan  (A*)  . 


tan  ( A1 ) 
tan ( A2) 
L1(XJ 
L2(X)  = 

RKX)  = 
R2(X) 


=  2  .  /  1  .  tan30 
=  2 ,  /  1 .  *  tan  120  = 

-  i  .  799  -  3 . 4641  »  (  X  4- 
1 . 799  +  1 . 1547  »  (X  + 

-  1 . 799  +  1 .  1547  »  (X  + 
1 .799  -  3.4641  #  (X  + 


1 . 15470 
-3 . 46410 
0.156 ) 
0.042) 

0 . 158) 

0 . 042) 


(2)  Ph  integral . 


Ph 


XI  L2(X) 

I  ! 

!  f (X)  !  f (Z)  dZ  dX 

I  I 

X3  Ll(X) 


XI 


X2  L2(X) 

I 

f  CX)  !  f (Z)  dZ  dX 
! 

R1  (X) 


X4  H2(X) 

I  I 

+  I  f (X)  I  f (Z)  dZ  dX 

I  I 

X2  R1 (X) 
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The  middle  double  integral  will  be  evaluated.  First  a 
substitution  U=  F(X) ,  dU  =  f (X)  dX,  is  made.  This  changes  the 
integral  to: 

U2  L2(X) 

I  f 

!  !  f(2)  dZ  dU 

!  I 

U1  Rl(X) 


and  changes  the  outer  integral  limits  to : 


U1 

F  (XI ) 

F(-0.158) 

43723 

U2  = 

F(X2) 

=  F(-0.042)  = 

48325 

Ph 

Ph 

be  low 

below 

X  =  2t  = 

Zb  = 

top 

bottom 

U  interval 

U 

F-KU)  L2(X) 

R1  (X) 

F(Zt) 

F(Zb) 

[  . 43723 .  . 44  ] 

.  435 

-0 . 16366  1 . 65852 

-1 . 60554 

.95139 

03550 

C  .  44  . 

.45] 

.  445 

-0 . 13830  1 . 68780 

-  1 . 77625 

.95428 

03785 

[  .  45  . 

.  46  ] 

.  455 

-0 . 1 1304  1 .71697 

-  1  . 74708 

.95701 

0403  1 

[  .  46  , 

.  47] 

.  465 

-0.08785  1.74606 

-  1  . 7  1700 

.95960 

04299 

L  .  47  , 

.  48  ] 

.  475 

-0.06271  1.77508 

-  1 . 68897 

.96206 

0456  1 

[  .  48  , 

.  48325  ] 

.  485 

-0.03761  1.80407 

-1  . 65999 

Q  /I  *7  Q 

0434G 

Ph  between 

?h  in 

top  and  bottom  Ph 

1  in 

rectangular 

U  interval 

F(2t)  -F(2b)  U 

interval 

interval 

[  .43723 ,  . 44] 

.9159 

00277 

.002537 

[  .  44  , 

.  45] 

.9164 

01 

.  009164 

[  .  45  , 

.  46  ] 

.9167 

0  1 

. 009167 

[  .  46  . 

.  47] 

.9166 

01 

. 009166 

[  .  47  , 

.  48  ] 

.9164 

0  1 

.  009164 

C  .  48  , 

.  48325  ] 

.9159 

00325 

. 002977 

Ph  =  .04218 

The  Ph  for  the  complete  target  is  given  in  paragraph  6b. 
d.  Angle  attack. 

( 1 )  Target  area . 

(a)  X-Z  plane . 

Hull,  Paragraph  6c(l)  (a)J^  has  the  computed  hull 
coordinates . 


2.  Turret. 


X6 

s 

M 

w/2 

S 

cos30» . 

25-sin30»t' .  75 

1 

Z6 

1/2 

sin30» . 

25  +  cos30^^ .  75 

5 

0. 

1585 

m 

0 

7745 

m 

X7 

M 

-w/2 

s 

cos30» ( 

- . 25) -Sin30 

»  . 

Z7 

1/2 

sin30* ( 

- . 25) +cos30 

. 

— 

- 

0 

5915 

m 

0 

5245 

m 

xa 

M 

w/2 

s 

cos30» . 

25-sin30^♦  (-. 

75) 

Z8 

-1/2 

Sin30» . 

25  +  cos30<»  (  -  . 

75) 

— 

0 

5915 

m 

- 

0 

5245 

m 

3 .  Barr i er . 


XIO 
Z  10 


M 


- (W/2^H)  I 
- (L/2  +  K)  I 


I  cos30’♦(-2.75)-s^n30^♦(-3.25) 
I  sino0s{-2.75)-^cos30»(-3.2t) 


-0.7567  m 
-4 . 1896  m 


(b)  Y-2  plana. 

1 .  Y  coordinates . 


Y1 

Y2 

Y6 

YIO 

YO 


-H 

-1.75  m 
0  . 
h 

2  m 

-H+Hb  =  -1.75+1.5 

-.25  m 

(h  +  YIO)  /2  =  (2.  -.25)  /2. 

0.875  m 


2_.  Z  coordinates.  See  paragraph  6d(l)  (a)  . 
(c)  X-W  plane . 

J_.  X  coordinates.  See  paragraph  6d(l)  (a)  . 
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W  coordinates.  W%  =  Z%  *  sin(  )  Y%  >  cos  (  ) 

W1  =  -1.799  »  .Sin30  -  1.75  *  cos30 

=  -2.4150  m 

W2  =  1.799  »  sinSO  +  0.  »  cos30 

=  0.8995  m 

W6  =  0.7745  #  sin30  +  2.  »  cos30 

=  2.1193  m 

WIO  =  "4.1896  *  sin30  -.25  *  cos30 

=  -2.3113  m 

WO  =  0.  *  Sin30  +  0.875  *  cosSO 

=  0 . 7578  tn 

Angles . 

D1  =  sqrt  (cos(Al)  +  (sin(  )  »sin(Al))  **2) 

=  sqrt  (cos30  +  (sin30  *  sin30)  »#2) 

=  0.90139 

cos(Bl)  =  cos(Al)  /D1  =  COS30  /O. 90139 

=  .9608 

sin(Bl)  =  sin(  )  »sin(Al)  /D1 

=  sin 30  *  sin30  /O. 90139 
=  0.27735 

tan(Bl)  =  sin<Bl)  /  cos(Bl) 

=  0.2837 

D2  =  sqrt  (cos(A2)  *»2  +  (3in(  )  «sin(A2)  ^^2) 

=  sqrt  (cosl20  **2  +  (sin30  <>slnl20)  »»2) 

=  0.66144 

cos(B2)  =  Cos(A2)  /D2  =  cosl20  /O. 66144 

=  -0.7559 

sin(B2)  =  sin(  )  »sin(A2)  /D2 

=  Sin30  »*sinl20  /O, 66144 
=  0,65465 

tan(B2)  =  sin(D2)  /  cos(B2) 

=  -0.36605 

Barrier  coverage. 


U-V  coordinates. 


M 

COS ( U ) 

sin (U) 

= 

sin (B2) 

-  c  o  s  ( B  2  )  1 

cos (V) 

3  in (V) 

-  s  i  n  ( B  1 ) 

cos ( B  1 )  1 

0 .65465 

0 . 75591 

-0 . 27735 

0 . 96081 

UlO 

=  M  1  X  1  (; 

-  .65465»(- 

.  7  567)  +  .7559t^(- 2. 3113)1 

vio 

1  WIO 

- . 27735»  (  - 

.7567)  +  .  9608it  (-2.3113)1 

=  I  -2 . 2425  m 

I -2 . 0 108  m 
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U1 

=  M 

XI 

sz 

. 65465* (-0 .  116)+. 75 59 *(-2.415) 

VI 

Wl 

- . 27735* (-0.116)+.9608*(-2.415) 

-1 .9014  m 
-2.2882  m 


b.  Coverage.  UlO  <  U1  and  VIO  >  VI.  Case  equals  3. 
Bight  side  of  barrier  covers  hull  and  intersects  with 
left  side  of  hull  at  point  11. 

Kll  =  (tan(B2)  »  XI  -  tan(Bl)  *  XIO  +  WIO  -  Wl) 

/  (tan(B2)  -  tan(Bl)l 

-  (-0.866  »(~0.116)  -  0.2887  *(-0.7567)  - 

2.3113  +  2.415)  /  (-0.866  -  0.2887) 

=  -0.3660  m 

Wll  =  (Wl  /  tan(B2)  -  WIO  /  tan(Bl)  +X10  -  XI) 

/  (1  /  tan(B2)  -  1  /  tan(Bl)) 

=  (-2.415  /(-0.866)  +  2.31)3  /O. 2887  -  0.7567  + 

0.116)  /  (1  /(-0.866)  -  1  /0.2887) 

=  -2. 1985  m 


5_.  Turret  coverage. 


a.  U-V  coordinates. 


U6 

V6 


M 


X6 

= 

W6 

.65465<K“0  .  1  585)  +  .  7559*2  .  1  2  93 
-  .  27735*t  (-0.  2  585) +.9608*2.  1  193 


1  .  4982m 
2 . 0802m 


U2 

=  M  1X2 

s 

.65465*0.  116+. 7559*0. 0995 

V2 

1  W2 

-.27735*0.  1  16+ . 9608*0 . 8995 

1  ^■\  rr  C  c  —  I 
\^  .  /  ^  ^  V  ill 

0 . 8321m 


b.  Extention  over  hull.  U6  >  U2  and  V6  >  V2 .  Case 

equals  2.  The  turret  area  may  extend  above  both  left 
and  right  edges  of  hull  area.  X7  <  X6  <  X8 : 
therefore,  the  turret  area  does  extend  above  both 
edges  of  hull  area.  Make  bottom  point  5  of  turret 
ax-ea  the  top  point  2  of  hull. 


(2)  Ph  integral.  The  Ph  is  given  in  paragraph  7b. 
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7.  Cornputer  ppo^ram  Top  ph .  Computer  program  'Top_Ph'  uses  the 
new  methodology  in  DFWP  to  compute  Ph  for  sample  weapons  against 
sample  targets.  Inputs  are  the  target  dimensions,  aspect  angle, 
attack  angle,  and  probability  distribution  data  for  the  weapon. 

Program  Top_Ph  uses  two  methods  for  the  integration,  one 
being  the  method  in  DFWP.  For  level  attack,  the  methods  are 
table  lookup,  the  method  in  DFWP,  and  the  numerical  method.  The 
numerical  method  in  DFWP  is  done  by  subroutine  SUMMATE .  For 
angle  attack,  the  integration  methods  are  done  by  subroutines 
SUMMATE  and  TRAPEZOID.  TRAPEZOID  does  the  integration  by  using 
the  trapezoid  method  and  without  the  change  of  variable. 
TRAPEZOID  is  made  very  accurate  by  doing  10,000  interv.al3. 


APPENDIX  A 
TOP  PH  OUTPUT 


HULL  WIDTH,  LENGTH,  and  HEIGHT? 

2,00  3.00  1.75 

TURRET  WIDTH.  LENGTH,  and  HEIGHT? 

0.50  1.50  2.00 

BARRIER  HEIGHT? 

1 . 50 

ASPECT  ANGLE  and  ATTACK  ANGLE? 

30  .  0  . 

XBIAS,  XDISPERSION? 

0.20  2.00 

WBIAS,  WDISP? 

0.00  1.00 

LOOKUP  SUMIvIATE 
PH  HULL  =  0.03500  0.3500 

PH  TURRET  =  0.15714  0.15714 

PH  =  0.19214  0.19214 
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HULL  WIDTH,  LENGTH,  and  HEIGHT? 


2.00  3.00  1.75 

TURRET  WIDTH,  LENGTH,  and  HEIGHT? 


0 . 50 

1 . 50 

2 . 00 

BARRIER  HEIGHT? 

1 . 50 

ASPECT 

ANGLE  and 

ATTACK  ANGLE? 

30 .  ! 

30  . 

XBIAS , 

XDISPERSION? 

0.20 

2 . 00 

WBIAS . 

WDISP? 

0 . 00 

1 . 00 

SUMMATE 

TRAPEZOID 

PHI  = 

0 . 14498 

0 . 14499 

PH2  = 

0 . 042  19 

0 . 04219 

PH.3  = 

0  .  15423 

0.15422 

PH  = 

0 . 34140 

0.34140 

HULL  WIDTH.  LENGTH,  and  HEIGHT? 

2 . 00 

3 . 00 

1 . 75 

TURRET 

WIDTH,  LENGTH,  and  HEIGHT 

0 . 50 

1 . 50 

2 . 00 

BARRIER  HEIGHT? 

1  .  50 

ASPECT 

ANGLE  and 

ATTACK  ANGLE? 

30  . 

30  . 

XBIAS . 

XDISPERSION? 

0 . 20 

2.00 

WBIAS , 

WDISP? 

0.00 

1 . 00 

SUMMATE 

TRAPEZOID 

HULL 

PHI  = 

0 . 1 1273 

0,11 273 

PH2  = 

0 . 02490 

0 . 02490 

PH.3  = 

0 , 08920 

0,08919 

FH 

0 . 22683 

0 . 22682 

TURRET 

PHI  = 

0.03332 

0 . 03332 

PH2  ~ 

0.01912 

0.01913 

PK3  = 

0 . 03248 

0 . 03248 

PH  = 

0 . 08492 

0 . 08493 

TARGET 

PH  = 

0.31175 

0.31  174 
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APPENDIX  B 


TOP  PH  CODE 


PROGBAM  PH 

REAL  Lh.  Lt 

DATA  DTR  /. 01745329/ 

COMMON  /TA/  XBOT,  WEOT ,  XTOP ,  WTOP ,  XMEAN,  XOISP, 
1  WDISP 

COMMON  /ROT/  C0S_A1 ,  SIN_A1 


C  ** 


I 


INPUT 

WRITE  (») 
WRITE  (i. 
READ  (*  ,’ 
WRITE  (1. 
WRITE  (») 
WRITE  (1. 
READ  ( » , ’ 
WRITE  (1. 
WRITE  (») 
WRITE  (1, 
READ  (#  ,  ’ 
WRITE  (1, 
WRITE  (*) 
WRITE  (1, 
READ  ( *  ,  ’ 
WRITE  (1, 
WRITE  (*) 
WRITE  (1, 
READ  (  »  ,  ’ 
WRITE  (1, 
WRITE  (») 
WRITE  (1, 
BEAD  (» , ’ 
WRITE  ( 1 , 


’  HULL  WIDTH,  LENGTH,  and  HEIGHT"^ 

’(A)’)  ’  HULL  WIDTH.  LENGTH,  and  HEIGHT? 

(3F7 . 0) ’ )  Wh ,  Lh ,  Hh 
’  (3F7 . 2)  ’ )  Wh .  Lh .  Hh 
’  TURRET  WIDTH,  LENGTH,  and  HEIGHT? 

’(A)’)  ’  TURRET  WIDTH,  LENGTH,  and  HEIGHT? 

(3F7 . OJ  ’ )  Wt .  Lt .  Ht 
' (3F7 . 2) ■ )  Wt .  Lt ,  Ht 
'  BARRIER  HEIGHT? 

’(A)’)  '  BARRIER  HEIGHT? 

(F7 . 0) ' )  Hb 
’  (F7 . 2)  ’  )  Hb 

’  ASPECT  ANGLE  and  ATTACK  ANGLE? 

'(A)’)  ’  ASPECT  ANGLE  and  ATTACK  ANGLE? 

(2F5 . OJ  ’ )  A1 .  TA 
’ (2F5 .0) ’ )  A1 ,  TA 
’  XBIAS,  XDISPERSION? 

'(A)’)  '  XBIAS,  XDISPERSION? 

(2F7 .0)  ’ )  XBIAS ,  XDISP 
' (2F7.2) ')  XBIAS.  XDISP 
’  WBIAS,  WDISP? 

’ (A) ’ )  '  WBIAS ,  WDISP? 

(2F7 . 0} ’ )  WBIAS .  WDISP 
'  r2F7 .2)  '  )  WBIAS  .  WDISP 


C  »»  ANGLES 


cos_ 

A1 

= 

cos 

(DTR 

Al) 

SIN 

A1 

SIN 

(DTR 

Al) 

TAN~ 

A1 

= 

TAN 

(DTR 

Al) 

COS 

A2 

= 

-SIN 

A1 

SIN 

A2 

cos 

A1 

TAN" 

A2 

SIN_ 

A2  / 

C0S_A2 

COS, 

TA 

COS 

(DTR 

tf 

TA) 

SIN 

TA 

SIN 

(DTK 

TA) 

C  COORDINATES 
C  AIMPOINT 


WMEAN . 


27 


c 


c  *  # 


c 


xo  = 

YO  = 
20  = 
HULL 
CALL 
Y1  = 
CALL 
Y2  = 

■  CALL 
CALL 
TURRET 
CALL 
Y6  = 
CALL 
CALL 
BARRIER 
CALL 
YIO 


0  . 

(Ht  -  Hh  +  Hb)  /  2. 

0. 

ROTATE  (-Wh/2..  -Lh/2.,  XI,  21) 

-Hh 

ROTATE  (Wh/2..  Lh/2..  X2 ,  22) 

0. 

ROTATE  (-Wh/2.,  Lh/2..  X3 ,  23) 
ROTATE  (Wh/2..  -Lh/2..  X4 ,  23) 

rotate  (Wt/2..  Lt/2..  X6.  26) 

Ht 

ROTATE  (-Wt/2.,  Lt/2.,  X7 ,  27) 
ROTATE  (Wt/2.,  -Lt/2.,  X8 ,  28) 

ROTATE  (-Wh/2.-Hh,  -Lh/2.-Hh,  XIO, 
=  -Hh  +  Hb 


I 


210) 


C  LEVEL  ATTACK  »»»»»###  it  if 

IF  (TA  .EQ.  0.)  THEM 

XMEAN  =  XO  +  XBIAS 
WMEAN  =  YO  +  WBIAS 
C  If*  HULL 

PH_HULL  = 

1  (DFN  ((X4-XMEAN)  /XDISP)  -  DFN  ( (X3- XMEAN)  /XDISP) ) 

2  (DFN  ((0. -WMEAN)  /WDISP)  -  DFN  ((YIO-WMEAN)  /WDISP) 


It 


WBOT  =  YIO 
WTOP  =  0 . 

CALL  SUMMATE  (X3 ,  X4 ,  0.,  0.,  PH_HULL_SUM) 
itif  TURRET 

PH_TURR  = 

1  (DFN  ((X8-XMEAN)  /XDISP)  -  DFN  ((X7-XMEAN)  /XDISP))  ^ 

2  (DFN  ((Y6-WMEAN)  /WDISP)  -  DFN  ((0. -WMEAN)  /WDISP)) 

WBOT  =  0 . 

WTOP  =  Y6 

CALL  SUMMATE  (X7 ,  X8 .  0..  0.,  PH  TURR_SUM) 

WRITE  (1 ,  ’  (A/3(A,2F9.4/) )  ’) 

1  '  LOOKUP  SUMMATE’ , 

2  ’  PH  HULL  =’.  PH_HULL,  PH_HULL_SUM, 

3  '  PH  TURRET  =’.  PH_TURR,  PH_TURr"sUM, 

4  'PH  = ’ ,  PH  HULL+PH  TURR," 

5  PH  HULL  SUM+PH  TUBR  SUM 


C  If  If  If  It  It  If  If  It  If  If  STRAIGHT  DOWN  ATTACK  if  it  #  if  if  if  it  it  if  it 
ELSE  IF  (TA  .EQ.  90.)  THEN 
XMEAN  =  XO  +  XBIAS 
WMEAN  =  20  +  WBIAS 
XBOT  =  XI 
WBOT  =  21 
XTOP  =  X2 
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n  'I*  if>  <o 


WTOP  =  Z2 

C  DIVIDE  HULL  into  THREE  SECTIONS 

IF  (XI  -LT.  X2)  THEN 
XL  =  XI 
XR  =  X2 

TAN_A_  =  TAM_A1 

ELSE 

XL  =  X2 
XR  =  XI 

TAN_A_  =  TAN_A2 
ENDIF 

C  PH  on  LEFT  SECTION 

CALL  SUMMATE  (X3.  XL.  TAN_A2 .  TAN_A1 .  PH1_SUM) 
CALL  TRAPEZOID  (X3 ,  XL,  TAN_A2 ,  TAN_A1 .  PH1_TRAP) 
C  »»  PH  on  MIDDLE  SECTION 

CALL  SUMMATE  (XL,  XR .  TAN_A_ .  TAN_A_ .  PH2_SUM) 
CALL  TRAPEZOID  (XL,  XR ,  TAN_A_ ,  TAnIa_ ,  PH2_TEAP) 
C  **  PH  on  RIGHT  SECTION 

CALL  SUMMATE  (XR.  X4 ,  TAH_A1 ,  TAN_A2 ,  PH3_SUM) 
CALL  TRAPEZOID  (XR,  X4 ,  TAN  A1 ,  TAN_A2 .  PH3_TRAP) 
WRITE  ( 1 .  '  (A/4 (A,2F10. 5/) )  '  ) ~ 


SUMMATE 

TRAPEZOID'  , 

PKl  =  ’ 

,  PH1_SUM. 

PHI  TRAP, 

PH2  =  ’ 

,  PH2_SUM, 

PH2_TRAP, 

PH3  =  ’ 

.  PH3,SUM, 

PH3_TRAP, 

PH  =  ’ 

DTT  V  Yf  t# 

,  X  11  i  _0  Uifv 

T  PH2_SUM  PH5 

PHl.TRAP  +  PH2_TRAP  +  PH3,TBAP 

C  »»»#»«#»»»  ANGLE  ATTACK 
ELSE 

C  W  COORDINATES 

WO  =  YO  COS_TA  +  ZO  »  SIN_TA 

W1  =  Y1  »  COS^TA  +  Z1  SInIta 

W2  =  Y2  *  COS’tA  +  Z2  »  SIN_TA 

W6  =  Y6  »  COS_TA  *  Z6  *  SIN^TA 

WIO  =  YIO  *  COS_TA  +  ZIO  »  SIN_TA 
C  »*  ANGLES 

D1  =  SORT  (COS_Al  »*2  +  (SIN„TA  »  SIN_A1)  **2) 

COS_Bl  =  COS..aT  /  D1 

SIN^Bl  =  SIN_TA  »  SrN_Al  /  D1 

TAN^Bl  =  SIN_B1  /  C0S_B1 

COS  V  =  -SIN_B1 

SIN_V  =  COS_Bl 

! 

D2  =  SORT  (COS_A2  **2  +  (SIN  TA  »  SIN  A2 '  **2) 

COS_B2  =  COS_A2  /  D2 
SIN_D2  =  SIN_TA  »  SIN_A2  /  D2 
TAV_,&2  =  SIN_B2  /  COS_B2 
cosju  =  SIN_i2 
SInIu  =  -COS_B2 
C  MEANS 

XMEAN  =  XO  +•  XBIAS 
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ooo  ooo  ooo 


WMEAN 


WO  +  WBIAS 


I 

C 

c 


c 


*»»»»  HULL 

**  POINT  1  and  BARRIER  POINT  in  U-V  COORDINATES. 
U1  =  XI  *  COS_U  +  W1  »  SIN_U 
VI  =  XI  »  COS  V  +  W\  SIN  V 


#  * 
*  * 


1 

1 


1 

1 


UlO  =  XIO  *  COS_U  +  WIO  »  SIN_U 

VIO  =  XIO  »  COS_V  +  WIO  »  SIN_V 
DETERMINE  SITUATION  of  BARRIER  COVERAGE  to  DETERMINE 
BOTTOM  POINT  of  EXPOSED  HULL. 

IF  (UlO  .GT.  U1  .AND.  VIO  . GT .  VI)  THEN 
BARRIER  PROJECTS  OVER  BOTH  BOTTOM  EDGES  OF  HULL. 

XBOT  =  XIO 

WBOT  =  WIO 

ELSE  IF  (VIO  .GT.  VI)  THEN 
**  BARRIER  PROJECTS  over  RIGHT  EDGE  of  HULL- 

COMPUTE  INTERSECTION  of  LEFT  EDGE  of  HULL  with  RIGHT 
**  EDGE  of  BARRIER. 

XBOT  =  (TAN  B2  *  XI  -  TAN_B1  *  XIO  +  WIO  •*  Wl) 

/  (TAN_B2  -  TAN_B1) 

WBOT  =  (Wl  y  TAN_B2  -  WIO  /  TAM_B1  +  XlO  -  XI) 

/  (l./TAN_B2  '  l./TAN_Bn 
ELSE  IF  (UlO  .GT.  Ul)  THEN 
»»  BARRIER  PROJECTS  over  LEFT  EDGE  of  HULL.  COMPUTE 
**  INTERSECTION  of  RIGHT  EDGE  of  HULL  with  LEFT  EDGE 
»»  of  BARRIER. 

XBOT  -  (TAN_B1  «  XI  -  TAN_B2  »  XIO  +  WIO  -  Wl) 

/  (TAN_Bl  -  TAN_B27 

WBOT  -  (Wl  /  TAN_B1  -  WIO  /  TAN_B2  +  WIO  -  XI) 

/  (l./TAN_Bl  -  l./TAN~B2) 

ELSE 

BARRIER  PROJECTS  under  HULL. 

XBOT  =  XI 


wB  0  T  =  w  1 
END  IF 
XTOP  =  X2 
WTOP  =  W2 


C  **  DIVIDE  EXPOSED  HULL  into  THREE  SECTIONS 

IF  (XI  .LT.  X2)  THEN 
XL  =  XI 
XR  =  X2 

TAN_B_  =  TAN_B1 

ELSE 

XL  =  X2 
XR  =  XI 

TAN_B_  =  TAN_B2 
ENDIF 

C  »»  PH  on  LEFT  SECTION  of  HULL. 

CALL  SUMMATE  (X3 ,  XL,  TAN_B2 ,  TAN_B1.  PH1_HULL_SUM) 
CALL  TRAPEZOID  (X3.  XL,  TAN_B2 ,  TAN_B1, 
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1  PHI  HULL_TRAP) 

»»  PH  on  MIDDLE  SECTION  of  HULL. 

CALL  SUMMATE  (XL,  XR .  TAN_E_ ,  TAN_B_ ,  PH2,HULL_SUM) 
CALL  TRAPEZOID  (XL,  XE ,  TAN_B_ ,  TAN_B_ , 

1  PH2_HULL_TRAP) 

**  PH  on  RIGHT  SECTION  of  HULL. 

CALL  SUMMATE  (XR,  X4 ,  TAN_B1.  TAN_B2 ,  PH3_HULL  SUM) 
CALL  TRAPEZOID  (XR,  X4 ,  TAN._B1,  TAN_B2, 

1  PH3_HULL_TRAP) 

**  PH  on  HULL 

PH_HULL_SUM  =  PH1_HULL_SUM  +  PH2^HULL_SUM  + 

1  PH3~HULL_SUM 

PH_HULL_TRAP  =  PHI  ._HULL_TRAP  +  PH2_HULL_TRAP  + 

1  PH3_HULL_TRAP  ' 

TURRET  ***** 

XBOT  =  X2 
WBOT  =  W2 
XTOP  =  X6 

WTOP  =  we 


**  POINTS  2  and 

6  in  U 

- 

V  COORDINATES 

U2  =  X2 

* 

COS_U 

+ 

W2 

SIN_U 

V2  =  X2 

* 

COS_V 

+ 

W2 

SIN,V 

U6  =  X6 

* 

COS^U 

+ 

W6 

SIN_U 

V6  =  X6 

* 

COS  V 

4- 

W6 

SIN^V 

TURRET  PROJECTS  over  BOTH  EDGES  of  HULL. 

IF  (U6  .GT.  U2  .AND.  V6  . GT .  V2)  THEN 
»»  DIVIDE  TURRET  above  HULL  into  THREE-  SECTIONS 


IF  (X2  .LT. 

X6)  THEN 

XL  =  X2 

XR  =  X6 

TAN_BOT 

=  TAN„B2 

TAN_TOP 

=  tan“bi 

ELSE 

XL  =  X6 

XR  =  X2 

TAN_B0T 

=  TAN_B1 

TAN~TOP 

•■=  TAN  B2 

END  IF 

PK  on  LEFT  SECTION  of  TURRET 
IF  (X7  .LT.  XL)  THEN 

CALL  SUMMATE  ( X7 , XL , TAN_B 1 , TAN_B 1 , 

1  PH1_TURR_SUM) 

CALL  TRAPEZOID  ( X7 . XL , TAN_3 1 , TAN  Bl, 
1  PHI  TURR  TRAP) 

ELSE 

PH1„TURR_SUM  =  0, 

PH1_TURR_SUM  =  0. 

ENDIF 

**  PH  on  MIDDLE  SECTION  of  TURRET 

CALL  SUMMATE  ( XL . XR , TAN_BOT , T AN_T0P , 

1  PH2  TURR  SUM) 
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1 


CALL  TRAPEZOID  ( XL , XR , TAN_BOT , TAN_TOP . 
PH2_TURR_TRAP) 

C  **  PH  on  RIGHT  SECTION  of  TURRET 

IF  (XR  .LT.  X8)  THEM 

CALL  SUMMATE  (XR.  X8 ,  TAN„B2 ,  TAN_B2 , 

1  PH3  TURR_SUM) 

CALL  TRAPEZOID  (XR,  XB ,  TAN_B2 ,  TAN_B2 , 
1  PH3_TURR_TRAP) 

ELSE 

PH3_TURR_SUM  =  0. 

PH3_TURR_TRAP  =0. 

EMDIF 


C  TURRET  PROJECTS  over  LEFT  EDGE  of  HULL. 

ELSE  IF  (V6  .GT.  V2)  THEN 

»»  COMPUTE  X  INTERSECTION  of  LEFT  EDGE  of  HULL  with 
**  RIGHT  EDGE  of  TURRET. 

XR  =  (TAN_B1  X2  -  TAN_B2  »  X6  +  W6  -  W2 ) 

1  /  (TAN_B1  -  TAN_B2) 

XR  =  MIN  (XR,  Xa) 

**  PH  on  LEFT  SECTION  of  TURRET 

CALL  SUMMATE  (X7.X6,TAN  Bl,TAN_Bl,  PH1_TURR  SUM) 
CALL  TRAPEZOID  (X7 . X6 . TAM_B 1 . TAM_B 1 , 

1  PHI  TURK  TRAP) 

»»  PH  on  RIGHT  SECTION  of  TURRET 

CALL  SUMMATE  (X6.XR,TAN  Bl.TAN  B2 ,  PH2  TURR  SUM) 
CALL  TRAPEZOID  (Xe , XR , TAN_B 1 , TAN_B2 , 

1  PH2,TURR_TRAP) 

n  * 

PH3_TURE_SUM  =  0. 

phs'turr'trap  =  0. 


C 


C 


C 


»»  TURRET  PROJECTS  over  RIGHT  EDGE  of  HULL. 
ELSE  IF  (U6  .GT.  U2)  THEN 


**  COMPUTE  X 

INTERSECTION  of 

R I GHT 

EDGE  of  HULL 

Tr  ^  w 

LEFT 

EDGE 

OF  TURRET 

XL 

=  (TAN_B1  » 

X6 

-  TAN 

B2  X2  +  W2 

-  we) 

1 

/' (TAN 

-  TAN_ 

B2) 

XL 

=  MAX  (XL, 

X7) 

PHI 

_TURR_SUM  = 

0  . 

PHI 

_TURB_TRAP 

=  0. 

PH  on  LEFT  SECTION  of  TURRET 

CALL  SUMMATE  (XL,  X6 ,  TAN_B2 ,  TAN  Bl, 

1  PH2_TURR_SUM) 

CALL  TRAPEZOID  (XL.  X6 .  TAN  B2 ,  TAN  Bl. 
1  PH2_TURE_TRAP) 

**  PH  on  RIGHT  SECTION  of  TURRET 

CALL  SUMMATE  (X6.  X3 .  TAM_B2 ,  TAM  B2 , 

I  PH3_TURE_SUM) 

CALL  TRAPEZOID  (X6,  X8 .  TAN_B2 .  TAN,B2 , 
1  PH3  TUBE  TRAP) 
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c  TUP.EET  IMAGE  IS  WITHIN  HULL  IMAGE 

ELSE 

PH1_TURR_SUM  =  0. 

phiIturr_trap  =  0. 

PH2,TURR_SUM  =  0. 

PH2_TURR_TP,AP  =  0. 

PH3_TURR_S’’M  =  0. 

PH3ItUER_TRAP  =  0. 

END  IF 

! 

PH_TURR_SUM  =  PH1_TURR_SUM  +  PH2_TURE_SUM  + 

1  ”  PH3_TURR_3UM 

PH_TURR_TRAP*=  PHl_TrjRR_TEAP  +  PH2_TURR_TRA?  + 

1  "  ph3Iturr_trap  T 

;VEITE  (1  ,  ’  (A/A/4(A.2F9.4/)  ,  A/ 4  ( A  ,  2F9 . 4/ )  ,A/A,2F9.4)  ’  ) 
1  '  SUMMATE  TRAPEZOID’, 

1  ’  HULL ’ , 

2  '  PHI  PH1_HULL  SUM,  PHI  HULL_TRAP , 

3  ’  PH2  ='.  PH2~HULL_’SUM.  PH2_HULL_TRAP , 

4  ’  PH3  =’.  PH3_HULL_SUM,  PH3..HULL_TRAP , 

5  ’PH  =’,  PH_HULL_SUM.  PH_HULL_TRAP , 

1  ’  TURRET’ , 

2  ’  PHI  =’,  PH1_TURR_SUM,  PH1_TURR_TRAP , 

3  ’  PH2  =’,  PH2'tURr7.SUM.  PH2_TURrItRAP  , 

4  ’  PH3  =’,  PH3^TURR_SUM,  PH3_TURR_TRAP , 

5  ’PH  =’,  PH_TUBR_SUM,  PH^TURR^TRAP , 

1  ’  TARGET’, 

2  'PH  =’ ,  PH  HULL  SUM  +  PH  TURR  SUM, 

3  PH_HULL_TRAP  +  PH_TURR_TRAP 
ENDIF 

STOP 

END 


SUBROUTINE  ROTATE  (X.  Z,  XX,  ZZ) 
common  /Rot/  cos_a,  sin_a 
XX  =  X  «  COS_A  -  Z  *  SIN_A 
ZZ  =  X  »  SIN_A  +  z  *  COS_A 
END 


FUNCTION  DFN  (X) 

1 

C  **  FROM  HASTINGS  APPROXIMATIONS  FOR  DIGITAL  COMPUTERS 

I 

F  =  0. 

AX  =-•  ABS  (X) 

IF  (AX  .GE.  5 . )  GO  TO  1 

F  =  (((((.  5383E-5<(AX+ .  488906E-4)  »AX+ .  380036E-4)  ^fAX  + 

1  .0032776263)  »AX+  .02  1  14  1006  1)  «fAX+.  049867  3469)  ^^AX+  1  . 

F  =  .5  /  (  (F^^^S)  »»2) 
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1  IF(X.GE.0,)F=1.  -F 

DFM  =  F 
END 


SUBEOUTINE  SUMMATE  (X_L.  X_R,  TAM_1,  TAiI_2  ,  PH) 


DIMENSION  DEVIATE(0 ; 99) 

COMMON  /TA/  X_1 ,W_1 .X_2 .W_2.  XMEAN , XDI SP , WMEAN ,  WPISP 
DATA  (DEVIATETk) ."k  ="o .  49)  / 


* 

-2 . 67583  , 

-2.  17010  . 

-1 . 95997 , 

-1 .81191 , 

- 1 . 69540 . 

-  1 . 59820  , 

-1.51411, 

-1 . 43954 . 

-  1 . 37220  , 

-1 . 31059 . 

* 

-  1 . 25356 . 

-  1 . 20036  , 

-1  .  15035 , 

-  1 .  10307 , 

-1 .05313 . 

-  1 . 01522 , 

-0 . 974 1 1 . 

-0 , 93459  , 

-0 . 89648 , 

-0 . 8596 1  , 

-0 . 82389 , 

-0 . 78919  , 

-0.75541 , 

-0 . 72248 . 

-0 . 69032 , 

-0 . 65884 . 

-0 . 62801  , 

-0.59776 , 

-0 . 56805 , 

-0 . 53884 , 

* 

-0.51008 , 

-0 . 48  173  . 

-0.45376 , 

-0 .42616 , 

-0 . 39S86  , 

* 

-0 . 37185 , 

~0 . 345  13  , 

-0.31864, 

-0 . 29238 . 

-0 . 26630  , 

-0 . 24043  , 

-0 . 21470  , 

-0 . 18913 . 

-0 .  16366  , 

-0 .  13830  . 

-0 .  1  1304  . 

-0 . 08785  . 

-0.06271 , 

-0 . 03761  , 

-0.01254  / 

DATA  (DEVIATE(K) ,  K  = 

50,  99)  / 

0 .01254  , 

0 . 0376  1  , 

0.06271 . 

0 . 08785  , 

0 .  11304  , 

0.  13830  , 

0 .  16366  , 

0.18913, 

0 . 2  1470  , 

0 . 24043  . 

* 

0 . 26630 . 

0 . 29238  , 

0 . 31864 , 

0.34513, 

0.37135, 

0.39886. 

0 . 42616  , 

0 . 45376  , 

0 . 48173  . 

0 ,51008  , 

« 

0.53884 , 

0 . 56805  , 

0 . 59776 . 

0.62801  , 

0 . 65884 , 

0 . 69032 , 

0.72248  , 

0 . 75541 . 

0.78919 , 

0 . 82389  . 

0 .8596 1 , 

0.89648. 

0.93459, 

0.97411 , 

1 .01522  , 

1 . 05813 , 

1 . 10307 . 

1 . 15035 . 

1 . 20036  . 

1 , 25356 ■ 

1 .31059  , 

1 .37220  , 

1 . 43954 , 

1.51411, 

1 , 59820  , 

1 . 69540  , 

1.81191, 

1 . 95997 , 

2.17010, 

2.57583  / 

C  *  *  Convert  coordinates  of  bottom  corner  Pi  and  top  corner  P2 
C  **  to  standard  deviations. 


XI  =  (X_ 

1 

-  XMEAN) 

/ 

XDISP 

W1  =  (W_ 

1 

-  WMEAN) 

/ 

WDISP 

X2  =  (X, 

2 

-  XMEAN) 

/ 

XDISP 

W2  =  (W 

2 

-  WMEAN) 

/ 

WDISP 

C  tn*  Convert  tangents  of  the  sides  to  standard  deviation  units. 
TANl  =  TAN_1  *  XDISP  /  WDISP 
TAN2  =  TAN_2  #  XDISP  /  WDISP 
C  Convert  left  boundary  to  standard  deviations. 

XL  =  (X_L  -  XMEAN)  /  XDISP 

C  Convert  to  probability  in  normal  distribution  table. 

UL  =  DFN  (XL) 

ULIOO  =  100.  ^  UL 

C  Convert  right  boundary  to  standard  deviations. 

XR  =  (X_R  -  XMEAN)  /  XDISP 

C  Convert  to  probability  in  normal  distribution  table. 

UR  =  DFN  (XR) 

URIOO  =  100.  *  UR 

C  **  I  is  subscript  for  inverse  normal  table  "DEVIATE'. 
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C  DEVIATE(O)  =  Standard  deviations  for  .005  probability 
C  DEVIATE(1}  =  Standard  deviations  for  .015  probability 
C 

C  DEVIATE(99)  =  Standard  deviations  for  .995  probability 
IF  (ULIOO  .NE.  100.)  THEN 
C  Compute  deviate  nearest  left  boundary. 

IL  =  ULIOO 

X  =  DEVIATE(IL) 

C  »»  WB  and  WT  are  bottom  and  top  of  rectangle  given  X. 

WB  =  W1  +  TANl  »  (X  -  XI) 

WT  =  W2  TAM2  *  (X  -  X2) 

**  Compute  subscript  for  1st  probability  in  'DEVIATE" 

»»  right  of  left  boundary. 

I  =  ULIOO  +  .5 

**  Initialize  summation  by  adding  or  subtracting  fractional 
part  of  beginning  increment. 

SUM  =  (DFN  (WT)  -  DFN  (WB) )  «  (I  -  ULIOO) 

ELSE 

SUM  =  0. 

END  IF 

C  »»  Compute  subscript  for  last  deviate  left  of  right  boundary. 
IMAX  =  INT  (URIOO  +  .5)  -  1 
DO  WHILE  (I  .LE.  IMAX) 

X  =  DEVIATE(I) 

WB  =  W1  +  TANl  ».  (X  -  XI) 

WT  =  W2  +  TAN2  *  (X  -  X2) 

SUM  =  DFN  (WT)  -  DFN  (WB)  +  SUM 
1  =  1  +  1 
END  DO 

IF  (URIOO  .NE.  100.)  THEN 
C  *»  Compute  deviate  nearest  right  boundary. 

IE  =  URIOO 
X  =  DEVIATE (IR) 

C  *+  Add  or  subtract  fractional  part  of  last  interval. 

WB  =  W1  +  TANl  (X  -  XI) 

WT  =  W2  +  TAN2  #  (X  •'  X2) 

SUM  =  (DFN  (WT)  -DFN  (WB)  )  »  (URIOO  -(IMAX+D)  +  SUM 

END  IF 

C  Multiply  by  interval  width  [1/100]. 

PH  =  SUM  /  100. 

END 


SUBROUTINE  TRAPEZOID  (X_L.  X_R .  TAN^l,  TAN_2 .  PH) 

I 

REAL  It  8  SUM 

DATA  SQRT_2PI  72.506628/ 

COMMON  /Ta7  X, 1 ,W_1 .X_2 . W_2 .  XMEAN , XD I SP , WMEAN , WD I SP 

I 

C  **  Convert  coordinates  of  bottom  corner  PI  and  top  corner  P2 
C  to  standard  deviations. 

XI  =  (X  1  -  XMEAN)  /  XDISP 
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W1  = 

(W  1  -  WMEAN)  / 

WDISP 

X2  = 

(X_2  -  XMEAN)  / 

XDISP 

W2  = 

(W~2  -  WMEAN)  / 

WDISP 

c 

Convert 

tangents  of  the 

sides  to  standard  deviation  units 

TAMl 

=  TAN_1  »  XDISP 

/  WDISP 

TAN2 

=  TAN_2  If  XDISP 

/  WDISP 

c 

Convert 

left  and  right 

boundaries  to  standard  deviations. 

XL  = 

(X_L  -  XMEAN)  / 

XDISP 

XR  = 

(X_R  -  XMEAN)  / 

XDISP 

c 

Compute 

interval  width 

for  10,000  intervals. 

H  =  (XE  -  XL)  /  10000. 

C  **  Start  trapezoid  summation. 

C  WB  and  WT  are  bottom  and  top  of  rectangle  given  X. 

WB  =  W1  +  TAMl  #  (XL  -  XI) 

WT  =  W2  +  TAW2  »  (XL  -  X2) 

SUM  =  EXP  (-XL»tXL/2.)  (DFM  (WT)  “DFN  (WB))  /  2. 
DO  T  =  1 . ,  9999 . 


X  =  T  #  H  -t-  XL 

WB  =  W1  +  TANl  »  tX  -  XI) 

WT  =  W2  +  TAN2  *  (X  -  X2) 

SUM  =  EXP  (-X»X/2.)  *  (DFN  (WT)  -DFN  (WB))  +  SUM 

END  DO 

WB  =  W1  +  TAMl  (Xfi  -  XI) 

WT  =  W2  +  TAN2  #  (XR  -  X2) 

SUM  =  EXP  ( ~XRi^XR/2 .  )  (DFN  (WT)  -  DFN  (WE))  /  2.  +  SUM 

PH  =  H  »  SUM  /  SQRT_2PI 

END 
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